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The Casimir friction problem for a pair of dielectric particles in 
relative motion is analyzed, utilizing a microscopic model in which we 
start from statistical mechanics for harmonically oscillating particles 
at finite temperature moving nonrelativistically with constant veloc- 
ity. The use of statistical mechanics in this context has in our opinion 
some definite advantages, in comparison with the more conventional 
quantum electro dynamic description of media that involves the use of 
a refractive index. The statistical-mechanical description is physical 
and direct, and the oscillator model, in spite of its simplicity, is nev- 
ertheless able to elucidate the essentials of the Casimir friction. As is 
known, there are diverging opinions about this kind of friction in the 
literature. Our treatment elaborates upon, and extends, an earlier 
theory presented by us back in 1992. There we found a finite fric- 
tion force at any finite temperature, whereas at zero temperature the 
model led to a zero force. As an additional development in the present 
paper we evaluate the energy dissipation making use of an exponen- 
tial cutoff truncating the relative motion of the oscillators. For the 
dissipation we also establish a general expression that is not limited 
to the simple oscillator model. 
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1 Introduction 



The Casimir effect - the original paper being Ref. [T] - has attracted a 
formidable increase of interest in the later years. For instance, as shown 
by a graph in the 2005 Lamoreaux paper [2], the number of citations per 
year to the 1948 Casimir paper has increased from about 10 to well above 
100 in the period from 1980 to 2010. 

The standard Casimir setup is that of two parallel metallic or dielectric 
plates at micron or sub-micron separations, the attractive force between them 
being calculated or measured. Later years have seen a considerable progress 
in the analysis of Casimir forces also between bodies of a more general shape, 
under static conditions. Some recent general treatments can be found in 
Refs. [31 SI El [6]. 

A special sub-class of Casimir phenomena is that of friction. The magni- 
tude of the Casimir friction force between plane surfaces in relative parallel 
motion has attracted considerable interest recently. The problem as such is 
not new; it was studied a long time ago by Teodorovich [7J and by Levitov 
[8]. Polarization currents fluctuating in the two bodies interact via the elec- 
tromagnetic field, transferring photons between the slabs. Because of the 
relative motion the photons are subject to Doppler shifts, and one should 
expect there to be a resulting dissipation of energy, meaning a friction force. 
The works mentioned were based upon macroscopic theory of the electromag- 
netic field in a medium, involving use of a refractive index. Most of the more 
recent literature has been formulated along the same macroscopic lines. Some 
of the later papers can be found in Refs. P [IUI [III [121 [131 [111 [ISl [IHl [III [II] . 

A complicating factor in the macroscopic theory is that there exists no 
natural rest inertial frame wherein the system of bodies is collectively at 
rest. This may be a major reason why diverging results are found in the 
current literature. For instance, Philbin and Leonhardt find there to be no 
quantum friction at all [13) , whereas most other papers find the friction force 
to be nonvanishing. Recent papers of Dedkov and Kyasov [17] and of Pendry 
[TB"] give useful overviews and comparison with earlier results. The title of 
Pendry's paper [IB] is illustrative for this kind of research: " Quantum friction 
- fact or fiction?" . 

Faced with this rather complex situation it becomes natural to inquire 
to what extent light can be shed on the problem by following an alternative 
approach: Instead of starting from quantum electrodynamics in a medium 
one can start with a simpler microscopic model in which only moving har- 
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monic oscillators are involved. Such an investigation is the main theme of 
the present paper. We shall consider two oscillators in relative motion hav- 
ing a constant nonrelativistic velocity. They are taken to represent a pair 
of polarizable particles. As we shall see, this oscillator model is in fact able 
to describe the essentials of the Casimir friction force. Some years ago we 
investigated such a microscopic model, for instantaneous [TS], as well as for 
non- instantaneous [20J, interactions. In view of the mentioned diverging re- 
sults obtained in this area of research, we find reason to consider this simple 
model anew, and to elaborate on it, the more so since results obtained in the 
model in essence should be valid for dielectric plates in relative motion also. 
The statistical mechanical approach implies that the Kubo formula [2TJ [22] 
will play an important role. This way of approaching the Casimir friction 
problem is apparently not so well known in the Casimir community, but we 
find it right to emphasize the microscopic method's usefulness. It leads to 
physical results without use of a heavy mathematical formalism. The micro- 
scopic approach has been followed in a recent investigation by Barton also 
(personal communication) . 

We assume in the following thermal equilibrium conditions, and assume 
the interaction between the oscillators to be weak. 

The results of the present paper can be summarized as follows: 

• There is a finite Casimir force, at any finite temperature T. 

• For the simple oscillator model the Casimir force vanishes at T — 0. 

• The Casimir energy dissipation is calculated, identifying the dissipation 
with the work done. Assuming infinite motion (i.e., the same velocity for all 
times), it is not clear how to distinguish between reversible energy change 
and irreversible energy dissipation. We handle this problem by introducing 
a convergence factor by means of which the interaction is limited to a finite 
time interval. Therewith the energy dissipation is calculated unambiguously. 

• The final expression for the energy dissipation, Eq. ( )27|) below, has a 
validity beyond the limitations of the harmonic oscillator model. Since this 
is a new result, it would be of interest to have it verified also in other ways, 
if possible. 

We mention finally that techniques similar to those used in the present 
letter have been used, for instance, in studies of the dynamical Casimir force 
associated with longitudinal motion of the plates (thus different from the 
lateral motion considered here). More information can be found in review 



3 



articles [2"3"| [21] , with further references therein. 

2 Calculation of the friction force 

We will consider the quantum mechanical two-oscillator system whose refer- 
ence state is the one of uncoupled motion corresponding to a Hamiltonian 
Hq. The equilibrium situation becomes perturbed by a time dependent term 
which we will write in general form as —Aq(t), where A is a time inde- 
pendent operator and q(t) a classical function of time whose explicit form 
depends on the specific properties of the system. The Hamiltonian becomes 
H = H — Aq{t). Moreover, we put 

- Aq(t) = ^(r(t)) Xl x 2 , (1) 

where -0(r) is the so-called coupling strength (i.e. ip is the classical potential 
between the oscillators). The separation between the oscillators is r, and 
X\ , x 2 are the internal vibrational coordinates of the oscillators. When the 
oscillators move with respect to each other the coupling has to vary in time. 
With nonrelativistic constant relative velocity v the interaction will vary as 

- Aq(t) = [V>(r ) + V^(r ) ■ vt + ...] Xl x 2 , (2) 

when expanded around the initial position r = ro at t = 0. The force between 
the oscillators, called B, is 

B = -(W(r))x 1 x 2 . (3) 

We ought here to mention the following point. In a mathematical sense the 
expansion (j2J requires vt to be small. Physically, we assume nevertheless 
Eq. to hold for all times, so that the interaction energy is taken to be 
proportional to t for all values of t. The natural opportunity of choosing 
q(t) = t in the interaction (pQ) thus has to be modified: as will be shown 
below, a convergence factor will be needed. 

Another point worth noticing is that the expression corresponds to 
first quantization only. Quantum electrodynamic processes such as emission 
and absorption of photons (second quantization) are not accounted for by 
the present model. They were considered, however, in Ref. [20] . 

The equilibrium situation with both oscillators at rest is represented by 
the first term in ([2]). It gives rise to a (reversible) equilibrium force. Thus 
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the friction must be connected with the second term. To simplify, we will for 
the moment neglect the first term. By this the two oscillators will be fully 
uncorrelated in their relative position r = ro- The friction force, due to the 
time dependence of the interaction (j2J), will be a small perturbation upon 
the equilibrium situation. This interaction leads to a response A(B(£)) in 
the thermal average of B. And this is where the Kubo formula, mentioned 
above, comes in [2T| 1221 125] : 

A(B(t)>= f cf) BA (t-t')q(t')dt' (4) 

J —oo 

(note that 4>ba means a vector), where the response function is given by 

fl A(*) = ^Tr{p[A,B(t)]}. (5) 

Here p is the density matrix and B(t) is the Heisenberg operator B(t) = 
e itH/h-Q e -ttH/n ^ w h ere b like A is time independent. Now with ([2]) and 
and with q(t) = t, expression (jSJ can be rewritten as 



with 



<P BA (t) = G0(£), (6) 

G = (V^)(v.VV), 
0(t) = Tr{pC(t)}, 

C{t) = —[x 1 x 2 ,x 1 {t)x 2 {t)]. (7) 
Thus with Eq. (j3J) and ^(t') = t' the force can be written as 

F = A(B(») = G f (j)(t - t')t'dt' = F r - 



roo 

Gt / <f>(u)du (8) 



where 



is part of the reversible force by which the part of the force that represents 
friction is 

F f = -G (f)(u)udu. (9) 

J 

Here the new variable u = t — t' has been introduced. The F r can be 
interpreted as a reversible force since it depends only upon position. This 
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interpretation is consistent with the result obtained for the dissipation in Sec. 
3 below; the F r will not contribute to the net total dissipation. 
If one again includes ^( r o) one has 

Gt = (W>)(v ■ Vtp)t (V^(ro) + tv ■ Vip + ...] 

= (V^(r + vi), (10) 

where r = ro+vt is the position at time t. By contrast, expression (j9]) changes 
sign when the velocity v changes sign, and it thus represents a friction force. 
(Observe that the velocity in (fIU|) merely represents the shift in position.) 
Equation is the same as result (2.11) in Ref. [H] 1, and the Fourier 
transformed version of it and Eq. (Q above as well, is 



du 



(11 



where 4>(u) = J °° (f)(t)e- iu)t dt (with <j>{t) = for t < 0). 

In Ref. [IH] the Fourier transformed version ffTT]) was used to obtain the 
explicit expression for the friction force. Here, we instead will use a different 
approach based on the expression (Q. As in the reference mentioned we 
then need the commutator (171). This entity again follows from the properties 
of quantized harmonic oscillators. We introduce annihilation and creation 
operators a and a) with commutation relations [a*, aj] = 1 (i — 1,2; other 
commutators vanish). As usual, cij(t) = aje~ luljt and a/(t) = a/e lWj *. With 
this the coordinates are 

x -=y^T^ + ^ < i2) 

where uii and U{ ii = 1,2) are the mass and eigenfrequency of each oscil- 
lator. To obtain <p{t) from (JTj) we first have to calculate {n\ri2\C(€)\Ti\n2), 
where \ri\n2) = \ni)\ri2) represents eigenstates with oscillators excited to lev- 
els n\ and n 2 . One has, when taking into account standard properties of the 
annihilation and creation operators, 

Li = (jii\a\ai(t) + aia\(t)\rii) = (2n« + 1) cos(u,it) +ism(uit), (13) 



3 Therc is a missing minus sign in Eq. (2.11) and the F r was not taken into account. 
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from which the thermal average follows after some computations as 

c(>(t) = ((r^lC^ln^)) 

= — (L\Li2 — -/vi L% 

in 2m,\ij}\ 2m 2 u; 2 

= D [(2(m) + 1) cos(wit) sin(u; 2 t) + (2(n 2 ) + 1) cos(w 2 *) sin^it)] . (14) 
Here 

D h 



With energy levels e n = (n + \)Tiuj the thermal average for the occupation 
numbers is 

2{m) + 1 = coth(~/3W). (15) 

In Ref. [19] expression ( 1T4|) was Fourier transformed to obtain the friction 
force as given by (JTTT) . As an alternative we will here use expression ( 1T4|) 
directly in Eq. fl9]). Then we get the integral 



f 

Jo 



te v cos(o>it) sin(o;2^)^ 



' A — ^(fc), *-o. (is) 



(») 2 + n?) 2 (i; 2 + si?,) 2 2Si 2 

Here f2i = o>i + u; 2 and f2 2 = oj\ — w 2 . As mentioned above a convergence 
factor e _r?i is needed, and the limit 77 — >- is taken. Then the f2 2 — term 
becomes a delta function with prefactor determined by the integral 

^ ;dx = l (17) 



oo(?7 2 + a: 2 ) 2 2 

From (I14p we also get this integral with u\ and oj 2 interchanged. This will 
then give the result ffT6l) with opposite sign with respect to the f2 2 -term. 
Adding up we get the difference of the prefactors 

1 ,1 swh.(\BKl 2 ) 

coth(-/3^i) - coth(-/3frw 2 ) 



'2^ ; v 2^ ; sinh(f/mcJi) sinh(|/3ftcc; 2 ) 
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sinh(i/3^a; 1 ) sinh(|/3^w 2 ) ' 



Multiplying ([To]) with ([18]) and including the factors D and G the friction 
force becomes 

which is also the result (3.14) of Ref. [19]. Again one notes that there is 
friction only when the oscillators have the same frequency, and ft should be 
finite, i.e. T > 0. 

When j3 — > oo, the expression (fT9l vanishes. According to the present 
oscillator model there is thus no friction force at zero temperature. An ob- 
jection against this result may be that it is somewhat singular due to the 
presence of the (5-function. Thus its physical significance may not be obvi- 
ous. However, r\ can be kept finite. This will smooth out the ^-function, and 
the f^-term in Eq. ( TT5"]) will give a contribution too. Note that this will not 
change our conclusions about a finite friction force for T > 0. (For finite rj, 
i.e. interaction like a short pulse, there will also be a contribution for T = 
due to the fii-term in Eq. (Ti~6]) .) But here we will assume r\ small by which 
the T = contribution can be disregarded. 

In Ref. [IS] the result ( I19p for the friction force was derived also by two 
other methods. These methods utilized the path integral formalism of quan- 
tum systems at thermal equilibrium [26] . The path integral can be identified 
with a classical polymer problem where imaginary time is a fourth dimension 
of length p. Thus the polymers stretch out in the fourth dimension and form 
closed loops of periodicity /3. For harmonic oscillators the correlation func- 
tion along the polymers is obtained in a straightforward way. With one of 
the methods the convolution of the correlation functions of both oscillators 
were needed. The resulting Fourier transform of this convolution was then 
identified with the response function <p(u) used in expression (TTTT) [25] . 

With the other method full thermal equilibrium was utilized. Then the 
relative motion of the oscillators was regarded as a harmonic oscillator motion 
with low frequency — > 0. Again with the path integral one can obtain 
the Fourier transform of the response function for the relative motion. The 
damping of the relative motion, that can be related to this response function, 
gives the friction force, and again the result (fT9|) was recovered. Thus the 
three methods used in Ref. [TJJ], as well as the modification considered in 
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the present paper, all lead to the same result, in contradiction to some other 
results in the literature, for instance that of Ref. [T3] . 

The result for the friction force was also extended to the situation with 
time-dependent or non-instantaneous interaction [20J. Then the full thermal 
equilibrium method was applicable to generalize the result. With the latter 
interaction there was also a friction from the self-interaction of a moving 
oscillator with itself. 

3 Dissipation of energy 

The presence of friction means dissipation of energy; the thermal energy of 
the system has to increase. However, this increase will be of second order in 
the velocity. As the Kubo formalism used above is limited to linear response, 
it is not immediately obvious how to evaluate the increase in energy or change 
in the Hamiltonian H . Simply replacing the operator B above with H will not 
work. However, for the present model it is possible to identify the dissipation 
directly with the work done. Work per unit time is force times velocity. 
During the perturbation period it is not obvious or possible how to distinguish 
between reversible change of energy and irreversible dissipation of energy. 
But this problem can be circumvented by considering the total energy change 
due to a perturbation that lasts for a finite time interval. This is the method 
that we will use in the following. 

Thus assume that the relative motion is finite in time, and that it starts 
at t = with maximum velocity v when the position is r = r . As t — > oo 
the motion is required to die out. To accomplish this we introduce the 
convergence factor e~ vt (77 — > 0) already used in Eq. ffl6l) . In the interaction 
(J2J), t is thus to be replaced with q(t) = te~ vt . As mentioned we consider 
only the time interval < t < 00, in which the velocity decays exponentially, 

v -> vi(t) = vq(t) = v(l - rfye^. (20) 

For rjt > , Vi(t) will now replace v in expression (fl9j) for the friction force. 
Altogether, the total energy dissipated will be 

/OO /"OO 
v 1 (t)-F f q(t)dt = vF f [q(t)] 2 dt=-v-F f , (21) 
-00 Jo 477 

where Ff is given by Eq. f fT9|) . Note that the reversible part of the force 
F r oc t — > q(t) as given by Eq. (jHJ) will not contribute to the dissipation since 
JS°q{t)q{t)dt = 0. 



9 



The above result for dissipation may be extended to more general cases 
for which q(t) shows an arbitrary variation with time. Thus q(t) need no 
longer be linked to a varying position in space as shown in Eq. ([I]). Under 
quite general conditions we may reinterpret the situation as one for which 
q{t) is associated with a "position", 

x{t) = q{t). (22) 

With this the "velocity" becomes 

v(t) = q(t). (23) 

As x(t) is a scalar quantity the corresponding operator B for the "force" will 
be 

B = ~{-Aq{t)) = A. (24) 
With relation (BJ the resulting "force" due to the perturbation becomes 

F f = f <p AA (t-t')q(t')dt\ (25) 

J — oo 

where 

<f>AA(t) = ±;Tr{p[A t A(t)]}- (26) 
With the associations made above the total dissipated energy becomes 



/oo roo r rt 

v(t)F f dt= / q(t) / <t> AA {t-t')q{t')dt' 



dt. (27) 



The expression ( 121]) is consistent with this more general result. In the integral 
( )2~Tj) a reversible part of the force was taken out. This corresponds to writing 

q(t') = q(t)-q(t)(t-t') + ... (28) 

The first term will not contribute to A£^, and we get (q(t) = for t < 0) 

AE d = - <j) AA {u)udu [q(t)]dt+ -, (29) 
J o Jo 



which is the result (J2TJ). This is seen by use of Eqs. (J3J), (JoJ) and (J9J). Then 
with g(t) = t — > te~ vt one has A = v ■ B and thus 0^ = v • <p BA = v • G <fi, 
by which — / °° (j) AA udu = v • F/. 
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The expression ( 127)) is thus the energy dissipation due to a perturbation 
—Aq(t) of the Hamiltonian. It is a general result, not necessarily limited 
to the simple model given by Eq. fll]). Since we used an indirect argument 
implying the net work done on an equivalent system, it would be of interest 
to verify this result by an independent method. The recent approach of 
Barton (personal communication) is in this context of interest, as he makes 
use of quantum mechanical perturbation theory to the second order in the 
interaction strength. If it could be shown that the different approaches lead to 
the same results for force and dissipation, it would be physically instructive. 
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